Droplet coalescence and breakup: Numerical 
simulations using moving mesh interface tracking 
with dynamic mesh adaptation 

Shaoping Quan and Jing Lou 
Institute of High Performance Computing, Singapore 138632 

October 18, 2010 
Abstract 

Droplet (s) breakup and coalescence have been simulated by the 
finite volume/moving mesh interface tracking method (MMIT) with 
adaptive mesh refining and coarsening. In this method, the interface 
is zero-thickness and moves in a Lagrangian fashion. The breakup 
and the coalescence of the interfaces are simulated by the mesh sep- 
aration and mesh combination schemes. Three cases are displayed in 
this video, which include the relaxation and breakup of an initially ex- 
tended droplet, two identical droplets off-center collision, and a droplet 
pair in a shear flow. The mesh adaptation is capable of capturing large 
deformations and thin regions, and the smallest length resolved is three 
orders of magnitude smaller than the droplet radius. The fluid dynam- 
ics videos of the simulations are to be presented in the Gallery of Fluid 
Motion, 2010. 



1 Introduction 

The two videos for numerical simulations of droplet (s) breakup and coales- 
cence are Videol and Video2. 

In the link videos, the simulations are performed using a finite vol- 
ume staggered mesh method coupled with moving mesh interface tracking 
(MMIT) scheme to locate the interface [2] . The advantages of this method 
is: 1) The interface is tracked by the triangular surfaces of the interior tetra- 
hedral mesh elements, and thus zero-thickness; 2) The jumps of the fluids' 
properties and stresses across the interface are implemented directly with- 
out smoothing; 3) The mass of each phase is naturally conserved due to 



the lagrangian motion of the interface. To capture the large deformations 
and wide range of length scales, a number of mesh adaptation schemes are 
implemented, and these mesh adaptation schemes distribute the mesh in an 
optimized way to achieve computing efficiency. The topological changes in 
multiphase flows are computed by the mesh separation and mesh combina- 
tion schemes [I]. 

Three cases are shown in the videos. The first case is the simulation of 
the relaxation and breakup of an initially extended droplet. The initial drop 
has a shape with a cylinder in the center and two bulbous ends. The surface 
forces are not balanced and thus the droplet experience a contraction and 
for this case the droplet eventually breaks up to three satellited drops. The 
density and viscosity ratios are 10, and 0.1, respectively, and the Ohnesorge 
number is 0.037. The result was reported previously [3], however, for the 
current simulation, the smallest neck radius is 10 times smaller than the 
previous simulation. The second case is an off-center droplet pair collision, 
and due to the large initial momentum, the two droplets coalesce and finally 
break up into five satellite drops. This case has a Weber number of 32, a 
Reynolds number of 80, and both the density and viscosity ratios of 10. The 
detail of the simulation can be found in [I] . The third case is a droplet pair 
in a shear flow with matched densities of the two fluids, and a viscosity ratio 
of 1.4, a Reynolds number of 0.1, and a Capillary number of 0.14. Initially 
the droplets are 5Rq (Rq is the radius of the spherical drop) away from 
each other streamwisely, and have a separation of 0.5i?o in the direction 
perpendicular to the free stream direction. The capability of our schemes 
in dealing with topological changes and resolving multiple length scales has 
been demonstrated. 
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